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MINIMUM AV, THREE-IMPULSE TRANSFER ONTO A 

TRANS-MARS ASYMPTOTIC VELOCITY VECTOR 

By William C. Bean 
Manned Spacecraft Center 

SUMMARY 


For the numerical study that is presented, the primary objective was to obtain 
a set of minimum AV, three -impulse maneuvers from a specified circular Earth 
parking orbit onto a constrained asymptotic velocity vector characteristic of Mars. 
A computer program based on an accelerated gradient method for function minimization 
was used. Twenty cases were examined in which the four parametric values (0°, 2°, 
4°, and 6°) for the out-of-plane angle of the asymptotic velocity vector and the five 
parametric values (8, 12, 24, 48, and 72 hours) for the period of the elliptical first 
transfer conic were used. It is believed that the numerical study makes the following 
quantitative contributions. 

1. The study verified numerically the existence of three-impulse maneuvers 
which provide measurable improvement in 18 of 20 cases studied, even for realistically 
small launch windows, over corresponding one-impulse maneuvers. For the cases 
examined, the improvement is from -0.23 percent to +7.40 percent lower velocity 
penalty, that is, total AV. 

2. For the cases examined, the study obtains the optimum relative distribution 
(among the impulses) of energy -changing and plane- changing duties. Previous studies 
have obtained minimum AV, three-impulse trajectories from sets of maneuvers in 
which all of the plane changing occurs at the high-altitude middle impulse. Imposition 
of this characteristic, which is valid for the minimum AV, infinite-time maneuver, 
would have been unrealistic for the (finite-time) cases examined, for which the low- 
altitude third impulse performed as much as 29. 3 percent of the plane changing. 

3. The study determines that the velocity penalty for the small- launch- window 
three-impulse cases examined is only 0. 38 to 5. 2 percent higher than for correspond- 
ing coplanar cases. 

4. The analysis discerns that the three-impulse maneuvers involve a 0. 13 to 

5. 8 percent greater velocity penalty than the corresponding absolute minimum AV, 
time-open three-impulse maneuvers involve. 


INTRODUCTION 


Studies frequently have emphasized energy considerations in the determination of 
a launch window for transfer from an Earth orbit onto a trans-Mars trajectory. These 
studies usually noted that a parking orbit could be established with an appropriate in- 
clination containing the required asymptotic direction. Therefore, the supposition was 
made that launch would occur with an instantaneous burn at the perigee of an escape 
hyperbola, with the plane of the departure trajectory coincident with the parking-orbit 
plane. 


For a specific mission, both the required asymptotic direction and the parking- 
orbit inclination vary in inertial space with time passage. Therefore, it is realistic to 
anticipate that the required asymptotic velocity vector may not lie in the parking- orbit 
plane. Consequently, the need to regulate fuel expenditure according to the dual re- 
quirements of performing plane changes and imparting energy indicates a further re- 
duction in launch windows for Earth-to-Mars missions. 

With regard to examination of the dual requirements, certain one- and two-burn 
requirements were studied by Deerwester (ref. 1). The following concepts were among 
the conclusions. 

1. The actual departure window is much smaller than the departure window 
dictated by energy requirements alone. 

2. The relationship between the declination of the required departure asymptote 
and the inclination of the parking orbit profoundly affects the size of the launch window. 

3. A quantitative comparison of out-of-plane escape maneuvers from circular 
orbit and, in particular, selection of a set of maneuvers yielding a minimum velocity 
penalty for a particular asymptotic velocity vector suggest extensive use of a digital 
computer. 

More recently, optimum three- impulse transfer onto an escape asymptote has 
been studied (ref. 2). For such a transfer, several additional degrees of freedom are 
permitted, including consideration of elliptical parking orbits. 

Finally, an analysis presented by Edelbaum (ref. 3) indicated the likelihood that 
a class of three-impulse maneuvers from a circle to a prescribed asymptotic velocity 
vector exists, each maneuver consuming less fuel than the corresponding minimum 
AV, one- and two-impulse transfers. The assumption is made that hyperbolic perigee 
can occur at an altitude lower than the circular parking-orbit radius. Also, it is 
suggested that, although unrealizably long flight times result, the lower bound for 
velocity expenditure is attained by allowing the altitude at which the second impulse is 
applied to increase without limit and by allowing the altitude at which the third impulse 
is applied to decrease to zero with respect to the origin of attraction. 

The conclusions of previous studies led to a numerical study of 20 cases in which 
the duration of the first transfer coasting arc and the target asymptotic velocity vector 
were constrained. In each case, the same target hyperbolic energy and the same 
lower bound for the perigee altitude of the escape hyperbola were fixed. The pair of 
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parameters for the study were prescribed values for the period of the coasting arc 

between the first and second impulses and for the out-of-plane angle D of the asymp- 
totic velocity vector. 

For each case, the objectives were to determine the following facts. 

1. A minimum AV, three-impulse maneuver from a specified circular Earth 
parking orbit onto the appropriate target asymptotic velocity vector 

2. The decrease in velocity expenditure relative to the appropriate reference 
one-impulse transfer 

3. The increase in velocity expenditure relative to the corresponding minimum 
AV, three- impulse coplanar maneuver which imparts the desired energy 

4. The relative distribution, among the three impulses, of both energy-imparting 
and plane-changing responsibilities 

5. The excess fuel expenditure over the corresponding absolute minimum AV, 
time-open, three-impulse transfer 

In the numerical study, a 262-nautical-mile Earth parking orbit was assumed. 
The hyperbolic energy was fixed by the requirement that the magnitude of the asymp- 
totic velocity vector V^ be 12 000 ft/sec. Hyperbolic perigee altitude was constrained 

in each case to be no lower than 100 nautical miles above the surface of the Earth. The 
set of values 8, 12, 24, 48, and 72 hours were considered for the period T^ of the 

first coasting arc; and the set of values 0°, 2°, 4°, and 6° were considered for the 

out-of-plane angle D of the vector V^. 

Finding a minimum AV maneuver involved a search in a class of pairs of trans- 
fer trajectories which could be ellipses, parabolas, or hyperbolas. However, specifi- 
cation of a period T^ of motion for the first transfer coasting path guaranteed 

ellipticity of this trajectory. Also, the numerical study avoided the following 
presumptions. 

1. That a given velocity increment was devoted entirely either to adding energy 
or to imparting plane changes 

2. That the three velocity increments for a three-impulse maneuver were 
mutually parallel or antiparallel 

3. That impulses must occur only at perigee or apogee of transfer conics 

4. That insertion occurred at perigee of the escape hyperbola 

5. That insertion must occur at or above circular parking-orbit altitude 
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For each case, a minimum AV, three-impulse maneuver subject to the specified con- 
straints was found by using a modified program of an accelerated gradient method for 
function minimization (ref. 4). This report summarizes and discusses the results of 
such a numerical study, with relation to a class of possible Earth-to-Mars missions. 

The author wishes to express gratitude to Ivan L. Johnson of the Mission Planning 
and Analysis Division, NASA Manned Spacecraft Center, for continued guidance in the 
adaptation of the program of the accelerated gradient method and for assistance in 
theoretical development. Also, the author wishes to thank T. N. Edelbaum and 
H. J. Kelley of Analytical Mechanics Associates for verbal communications and 
assistance. 


SYMBOLS 

_ k. _ k. n" 1 ^ 

[A] orthogonal 3 x 3 matrix relating (i, j,k) to (I, J, K) 

a semimajor axis of the terminal hyperbola, Earth radii (e. r. ) 

a. semimajor axis of the elliptical coast between the first and second 

impulses, e. r. 

/ -A — k Xfp — i. _k 'P 

[B] orthogonal 2 x 2 matrix relating ( ip, jl to (i, j) 

b semiminor axis of the terminal hyperbola, e. r. 

D declination of V w , deg 

d . V - , (e. r. )^/hr 

e eccentricity of the terminal hyperbola 

F performance index (the sum of the magnitudes of the velocity increments for 

a given maneuver), e. r. /hr 

F* augmented performance index, equation (41) 

g column vector of the original set of constraints with the components g^, g^, 

g 3 , and g 4 

g’ column vector of amended set of constraints with the components g^, g ^ , 

g 3 , and g 4 

gj intermediate equality constraint on a semimajor axis of the first coast 
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gg, gg terminal equality constraints on V M declination and semimajor axis, 
a respectively 

gg' amended declination constraint 
g^ terminal inequality constraint on the radius of perigee 

h scalar angular momentum of the terminal hyperbola, (e. r. )^/hr 

I unit vector from the center of the Earth to the position of the vehicle at 

initial time 


i, j, k orthonormal right-handed set of vectors 


V j P 


J 

K 

L 

M 

m 

N 

P 

r 

r 

r 

P 

— V. 

r 

P 

RA 


unit orthogonal vectors in the plane of the hyperbola 

unit vector from the center of the hyperbola along the asymptote 
coparallel to 

unit vector in the direction of motion of the vehicle at initial time 
unit vector such that K = I x J 
y(yw - zv) - x(zw - xw) 
x(xv - yu) - z(yw - zv) 

slope of the hyperbolic asymptote with respect to the ^f,, system 
z(zu - xw) - y(xv - yu) 

semilatus rectum of the terminal hyperbola, e. r. 
scalar magnitude of r, e.r. 
terminal position vector (at insertion) 
scalar magnitude of r^, e. r. 

radius -of -perigee vector of the terminal hyperbola 
right ascension, deg 



T 1 period of the elliptical coast between the first and second impulses for a 
given three-impulse maneuver, hr 

u, v, w velocity components in the (I, J, K) coordinate system of the target 
hyperbola at r 

-A 

V scalar magnitude of V , e.r./hr 

P P 


velocity vector at perigee for the terminal hyperbola 


V asymptotic velocity vector of the terminal hyperbola, e. r./hr 


x, y, z position components of r in the (I, J, K) coordinate system 
a control vector with the components a ^,a^, ... ,a ^ 


2 control vector a which minimizes F (a) 

/3..,/3 0 change in generalized eccentric anomaly between AV 1 and AV 0 and 
between AVg and AV 3 , respectively, (e.r.) 

AV- column vector of the initial velocity increment with the components Au. , 

1 AVj, and AWj = ag, and respectively, e.r. /hr 

AV„ column vector of the intermediate velocity increment with the components 
1 Au 2 , Av 2 , and AWg = a^, a Q , and respectively, e.r./hr 

AVg column vector of the final velocity increment with the components Au 3 , 
Av 3 , and Aw 3 = a 15 a 2> and a g, respectively, e.r./hr 

9 polar angle referenced from a specific axis, deg 

e true anomaly of Q 1 (fig. 1) with respect to the perigee of the first transfer 

conic, deg 

6 „ true anomaly of Q„ (fig. 1) with respect to the perigee of the first transfer 
conic, deg 

e 2 true anomaly of Qg (fig. 1) with respect to the perigee of the second 
transfer conic, deg 
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e g’ true anomaly of Qg (fig. 1) with respect to the perigee of the second 
transfer conic, deg 

e g" true anomaly of Qg (fig. 1) with respect to the perigee of the terminal 
hyperbola, deg 

X column vector of constant Lagrange multipliers with the components A 1 , 

Ag, and X 

3 2 

Ijl gravitational constant of the Earth, (e. r. ) hr 

<p argument of perigee with respect to the terminal position vector, deg 

Subscript : 

p conditions at the hyperbolic perigee 

Superscripts : 

* desired value at convergence 

T matrix transpose 

-1 matrix inverse 


METHOD 


Geometric symmetry permitted simplification of the mathematical model 

for minimum AV orbit transfer from a circle to a vector by substitution of an 

equivalent problem, namely, that of finding a minimum AV transfer, with the first 

impulse at a fixed point on the circular parking orbit, to a fixed magnitude and 

out-of-plane angle with respect to the parking-orbit plane. To solve the new problem, 

the in-plane angle for the vector was regarded as a parameter to be optimized in 

the solution. Then, in the first problem, any desired prespecified value for the in- 
plane angle could be regarded as having been obtained simply by a rotation through the 
appropriate angle of the initial point on the circular orbit. 

The solution was expressed in a right-handed, orthonormal inertial coordinate 
system (I, J, K) based on the parking-orbit plane. In this system, I is a unit vector 
directed from the center of the Earth to the position of the vehicle at initial time, J 
is a unit vector in the direction of motion of the vehicle at initial time, and K is a 
unit vector given by the cross product K = I x J. Because there was no loss of 
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generality, it was assumed that the parking orbit was in the equatorial plane with I 
directed toward Aries, K directed toward the North Pole, and J completing the 

right-handed system (I, J, K). Subsequently, the out- of -plane angle of the V M vector 

was referred to as declination and the in-plane angle as right ascension. 

For each of the 20 computer runs generated in the numerical study, the transfer 
maneuver obtained was required to minimize the performance index 


F(«) 


V 1 


i. 



t> 

< 

H - 4 

4 - 

AV 2 

+ 

AV 3 


( 1 ) 


subject to the intermediate constraint gj * 0, the terminal equality constraints 
g 2 = g 3 = 0, and the terminal inequality constraint g^ s 0, where the column vector 
g was defined as 


g l 


a l a i* 

g 2 


sin D - (sin D)* 

g 3 


a - a* 

_ g 4_ 


* 

1 

1 


( 2 ) 


Fixing (sin D)* sufficed to fix the declination angle D*, because -90° s D* s 90° 
by definition. Furthermore, a^* was defined so that a corresponding desired value 

for the period T^* was fixed according to the equation 


T 


j* = 2v 



(3) 


The value a j* was selected so that * = 12 000 ft/sec, where 


* = 



(4) 


The value r * was fixed at 1. 0290365646 e. r. (at a 100-nautical-mile altitude above 
P 

the surface of the Earth). 
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In using the program of the accelerated gradient method, it was useful to express 
g as a vector function of intermediate and terminal position and velocity variables. 
The resulting formula for g 2 involved fractions; therefore, subsequent hand calcula- 
tions of certain partial derivatives used in the program were simplified by the sub- 
stitution of gg’ for gg, where 


g 2 ’ = -erhg 2 


(5) 


It was preferable to use a formula for 


r , namely 


r _ h 2 (l + e) 
P M 


( 6 ) 


which is valid for all conics, including parabolas, because near -parabolic conics oc- 
curred as intermediate iteratives in the optimization process. 

For cases in which D* = 0°, it was found that inclusion of the second constraint 
introduced a redundancy with associated difficulties in obtaining convergence. The 
redundancy occurred because the optimization program haturally selects an in-plane 
solution in the absence of contrary constraints. Accordingly, for such cases, the 
constraint g 2 ' = 0 was deleted, leaving 


a 


1 



g = 


a - a* 



(7) 


with g, = g« = 0 and g 0 = r * - r < 0. 

°1 °2 3 p p 

Finally, consider the second constraint, g 2 = 0 initially, where 


g 2 = sin D - (sin D)* 


( 8 ) 


In spherical coordinates, the velocity vector at infinity is 


V 


oo 


Vjd cos RA cos D + J sin RA cos D + K sin 


D) 


(9) 


9 


Therefore 


(to) 


sin D = 


By analysis, for a hyperbola 


Ho = 


(V^p)~ 


(l + m 2 ) 


But for all conics 


1/2 


m 


= b = - a (e 2 - I. L. = -( e 2 - l) 


1/2 


Therefore 


= ( 1 + m 2 ) 


1/2 


Equation (11) becomes 


i = 

oo 


^ ( 2 1 ) 1 ^ 

lV U - lj . t 


However, by definition 


l = 

OO 


f 

P 7* 


t 3 

'p V p 


( 10 ) 


( 11 ) 


(12) 


(13) 


(14) 


( 15 ) 
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Therefore, equation (14) becomes 



( 16 ) 


Using the fact that 



(17) 


converts equation (10) into 


sin D = 



(18) 


As an intermediate step toward obtaining the necessary conditions for the prob- 
lem, it was important to express the constraint gg as a function of the components 

with respect to the initial rectangular coordinate system of the position and velocity 
vectors immediately after application of the third velocity increment. In other words, 
it was useful to change the right member of equation (18) from a function of state vari- 
ables at hyperbolic perigee to a function of state variables immediately after insertion 
into the hyperbola. This change can be made by using the following equations. 


r = xl + y J + zK 

> (19) 

V = ul + vJ + wK J 
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The unit vectors i, j, and k are defined as follows, where i and j lie in the 
hyperbolic plane of motion. 


k = 



r 

r 


\ 




3 = k x i = 



( 20 ) 


where h = |r x v| is the scalar angular momentum. 

Using equation (20), the following coordinate transformation can be written in 
matrix form. 


i 


I 

i 

= [A] 

J 

lx 


K 


( 21 ) 


where [ A] is the 3x3 square matrix given by 


x 

r 


y 

r 


z 

r 


[A] = 


z(zu - xw) - y(xv - yu) xfov - yu) - z(yw - zv) 
rh rh 


y(yw - zv) - x(zu - xw) 
rh 


(yw - zv) 
h 


(zu - xw) 
h 


(xv - yu) 
h 


( 22 ) 


The orthogonality of [ A] 


implies that 




therefore 


I 


i 

J 

= [a] t 

T 

K_ 


k 


( 23 ) 


12 



Accordingly, equations (19) become, after simplification 


where 



d = r • V = xu+yv + zw 


( 24 ) 


(25) 


The symbols N, M, and L are defined as 

N = z(zu - xw) - y(xv - yu) ' 
M = x(xv - yu) - z(yw - zv) > 
L = y(yw - zv) - x(zw - xw) 


(26) 


To convert equation (18) into an expression involving only quantities related to the 
initial rectangular reference frame at the instant of hyperbolic insertion, a pair of 

-fe -*■ -v -X 

rectangular coordinate axis systems (i, j) and (i , j ) in the plane of the hyperbola 

P P 

was defined. Therefore, it is possible to define an argument-of -perigee angle <p in 
such a manner that the (L,f ) system is related to the (i, j) system by 

r r 



(27) 


where [ B] is the 2x2 square matrix 


[B] 


cos <p sin 4> 

-sin 0 cos 0 


(28) 
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Therefore 


i 

p 


COS 0 

sin 0 


"(xl + yJ + zK) " 
r 

L 3 pJ 


-sin 0 

cos 0 


(NI + M J + NK) 
rh 


( 29 ) 


that is 


1 

i— >• L 

•° j 


Ai 



(xl + yJ + zK)cos <p 


• (xl + y J + zK)sin 0 
r 


(NI + MJ + LK)sin 0 
+ rh 

(NI + MJ + LK)cos <P 
+ rh 


(30) 


Equations (30) are equivalent to the pair of equations 


-r 1 (x cos 0 N sin <A T (y cos 0 M sin (p\ . 

*p \ r + rh / + \ r + rh / 


( z cos 0 h sin 0) 
+ \ r + rh / 

3 P " (" 


K 


x sin 0 N cos 0\ T ( y sin 0 M cos 0) T 
+ — ~ rh ) 3 


+ - 


r rh 

z sin 0 L cos 0 


rh 


K 


(31) 


The k- components of i and j are 

lr r 


p _ ( z cos 


0\ ( L sin 0^ 

) + l rh / 


W p _ (z cos <p\ ( L sin 4>\ 
V r / + \ rh ) 




(32) 
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Therefore, equation (18) becomes 


sin D = 


1 z cos 0 L sin 0 ( 2 

e r + rh " ' e 



z sin 0 L cos 0\ 
r + rh / 


( 33 ) 


It is useful to express e, h, cos 0, sin 0, and L in terms of x, y, z, u, 
v, and w. According to the laws of orbital mechanics 


h = [r 2 V 2 - (r . v) 2 ] ^ 


(34) 


Similarly 


e = 


1 - 



(35) 


The value of L was determined in equation (26). The remaining steps are to express 
cos 0 and sin 0 in terms of x, y, z, u, v, and w. 

The polar equation for a conic is 



h 


2 


+ e cos (-0)] 


(36) 


Therefore 


COS 0 



(37) 


2 2 

Because cos 0 + sin 0 = 1 and because d > 0 on the upper half of a conic and 
d < 0 on the lower half of a conic 


sin 0 = - 



(38) 
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Accordingly, it is helpful to replace the constraint g 2 by g 2 ' = -erhgg. It is 

possible to consider g ’ = 0, because g' = 0 if and only if g 9 = 0. From 
equation (35) 

gg' = |erh(sin D)* + zh cos <p + L sin <p 


7 2 V /2 

ye - 1 J (-zh sin 0 + L cos 0)J| 


(39) 


In all cases, the control vector a was defined as 


a - [a v q? 2 , a y a 4 , « 5 , "g, a 7 > a g> a 9’ ° 10’ "ll) 

= ^Au 3 , Av 3 , Aw 3 , /3 2 , Au 2 , Av 2 , Aw 2 , /3 1? AUj, Av^ AwJ T (40) 


In equation (40), /3 and /3„ were the changes in generalized eccentric anomaly 
between AV^ and AV 2 and between AV 2 and AVg, respectively. 

The initial position on the 262-nautical-mile parking orbit in each case was 
selected on the positive X-axis. Prior to the first transfer maneuver, the character- 
istic speed in the circular parking orbit of radius 1. 07607599156 e. r. was 25 003 ft/sec. 

Verification of a relative minimum for F was obtained by comparing the values 
found for the components of g' and SF' /da with zero, where 


F' 


T 

F + a g* 


(41) 


Therefore 


9F’~| _ /9F» 9F' 9F* \ T 

,9a J" ^9a 1 ’9a 2 ’"-’9« 1 J 


(0, 0, . . . , 0) T 


(42) 


f 


is a necessary condition for a minimum of F subject to the g- constraints. 

It was not possible to resolve a velocity impulse into components such that one 
component did only plane changing and did not alter the energy state. Nevertheless, it 
was possible to resolve an impulse into two components, one normal to the current 
orbital plane and the other lying in the orbital plane . Because only the normal com- 
ponent of an impulse resulted in orbital plane change, that component was called the 
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*■ 


* 


I, 


t- 


"component of the impulse associated with plane changing. " The in-plane component 
was called the "component of the impulse associated with energy changing." The 
"amount of plane changing accomplished by a three-impulse maneuver" was the sum of 
the magnitudes of the components associated with plane changing over all the impulses 
for the maneuver. The "percent of plane changing accomplished by a certain impulse" 
was defined as the percentage equivalent of the ratio of the magnitude of the "component 
of this plane -changing impulse" with respect to "the amount of plane changing associ- 
ated with the entire maneuver." A similar definition was adopted for "percent of 
energy changing accomplished by a certain impulse. " 


Each of the 20 computer runs with the accelerated gradient method program as- 
sumed an impulse-coast-impulse-coast-impulse type of maneuver. For each case, 
initial position and velocity components on a circular Earth parking orbit were input 
to the program. First guesses for a were made to start the program converging to 
find the value of a* resulting in a minimum total velocity expenditure. The geometry 
for a typical minimum AV, three-impulse transfer sequence is shown in figure 1. The 


first impulse AV^ for the maneuver was 

made on the parking orbit at a true anomaly 
0 j beyond the perigee of the first transfer 

conic. The second impulse AVg was 

made at a true anomaly e ^ beyond the 

perigee of the first transfer conic, causing 
injection into a second transfer conic. The 
position of the second impulse occurred at 
a true anomaly e beyond the perigee of 

the second transfer conic. Injection into a 
hyperbolic escape trajectory then occurred 
at a true anomaly e ^ beyond the perigee 

of the second transfer conic, that is, at a 
true anomaly e ^ ' beyond the perigee of 

the target hyperbola. 

For each of the in-plane cases, 
where D* = 0, the minimum AV transfer 
sequence was found to be a "double 
Hohmann" type of maneuver, that is, a 
180° -180° perigee-apogee-perigee type of 
maneuver. Specifically, for each in-plane 
case, e ^ = e g* = 0 g" =0° and 0 g = 

0 2 * = 180°. Or, for each in-plane case, 



TRAJECTORY 

TRAJECTORY 

TRAJECTORY 


0 CIRCULAR EARTH PARKING ORBIT 

1 FIRST TRANSFER ELLIPSE 

2 SECOND TRANSFER ELLIPSE 


TRAJECTORY 
Q l' Q 2' °3 


3 TERMINAL HYPERBOLA 

POINTS OF APPLICATION OF FIRST, SECOND, AND 
THIRD VELOCITY INCREM£NTS 
FIRST, SECOND, AND THIRD VELOCITY 
INCREMENTS 

APOGEE POINTS OF FIRST AND SECOND TRANSFER 
CONICS 


p. , p_, P, PERIGEE POINTS OF FIRST TRANSFER CONIC, 

1 * J SECOND TRANSFER CONIC, AND TERMINAL 

HYPERBOLA 


Figure 1. - Geometry for a typical mini- 
mum AV, three -impulse transfer 
sequence. 


the first impulse was at the perigee of the first transfer conic, the second impulse was 
at the apogee of both the first and the second transfer conics, and the third impulse was 
at the perigee of both the second transfer conic and the terminal hyperbolic trajectory. 
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For all in -plane and out-of -plane cases, except the two cases in which 
Tj = 12 hours and D = 2° and in which = 8 hours and D = 2°, the inequality con- 
straint was satisfied as if it were an equality constraint; that is, the minimum permis- 
sible value for target hyperbolic perigee r * was invariably attained. For the two 

P 

exceptional cases, r exceeded r * by 0. 04697 and 0. 03731 e. r. , respectively. 

P P 

Furthermore, in all cases, a retrofire was accomplished at the intermediate impulse. 

(The vector AV 2 always had a negative component in the direction of motion. ) The 

Hohmann features of the minimum AV transfer maneuvers apply only to in-plane se- 
quences. Computer rim 4 (T^ = 8 hours and D = 6°) is a sufficiently typical transfer 

sequence for a full three-dimensional case and is illustrated in figure 1. 

In figure 1, trajectory 0 is the circular parking orbit about the Earth. Trajec- 
tories 1, 2, and 3 are the first and second transfer ellipses and the terminal hyper- 
bolic path, respectively. The actual transfer flight path is from to Q 2 on 

trajectory 1, from Q 2 to Qg on trajectory 2, and thereafter on trajectory 3. The points 

Qj, Q 2 , and Q g are the points of application of the first, second, and third velocity 

increments, respectively. The points P^, AP^, P 2 , and AP 2 are the perigee and 

apogee points of the first transfer ellipse and the perigee and apogee points of the sec- 
ond transfer ellipse, respectively. The point P^ is the perigee of the terminal hyper- 
bolic trajectory. Impulses do not occur at apsidal points. In fact, the points P^, AP 2 , 
and P 2 are not even reached in this illustrative minimum AV maneuver. 

RESULTS 


Illustrated in figure 2 is the fact that the three-impulse technique provides some 
improvement in launch window over the reference one-impulse method, but falls short 
of the absolute minimum AV maneuver (corresponding to infinite T^) referred to in 
reference 3. 

As an illustration of the advantage of the three-impulse technique, note that for 
D = 6°, the minimum AV, three-impulse maneuver for T^ = 72 hours requires a 

velocity expenditure of 2. 15203 e.r./hr, a 7.4 percent smaller velocity penalty than 
for the reference one-impulse case which requires a total velocity increment of 
2.32409 e.r./hr. 

The velocity expenditure reduction attainable by increasing the period T^ of the 

first coast ellipse for a given target declination D is comparable. For D = 6°, the 
fuel expenditure for the case T^ = 72 hours is 2. 15203 e. r./hr, only 3. 9 percent less 

than for the most expensive minimum AV case examined for D = 6°; namely, 

T j = 8 hours, for which the velocity expended is 2. 23849 e. r. /hr. 
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DECLINATION OF THE ASYMPTOTIC VELOCITY VECTOR. DEG 


Figure 2. - Minimum AV velocity ex- 
penditure as a function of D (with 
a parameter T^). 


For the cases examined, the non- 
coplanarity velocity penalty is measurable, 
but not drastic (fig. 3). To illustrate the 
magnitude of the plane-change penalty, con- 
sider Tj = 12 hours. The velocity expend- 
iture for the D = 6° out-of -plane case is 
is 2. 21599 e. r. /hr, or 4. 3 percent higher 
than for the corresponding D = 0° co- 
planar case, with a velocity expended of 
2. 12430 e. r. /hr. For T^ = 72 hours, the 

velocity requirements for the 6° case ex- 
ceed those for the 0° case by only 
1. 6 percent. 

Another result is that for the cases 
considered, the typical minimum AV, 
three-impulse maneuver does not come 
close to completely uncoupling the energy- 
changing and plane -changing duties among 
the three impulses (table I). As a particu- 
lar example, for computer run 4 
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Figure 3. - Velocity expenditure as a function of 

for a three -impulse sequence. 
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TABLE I. - DISTRIBUTION OF PLANE-CHANGING AND ENERGY-CHANGING 


DUTIES FOR THE MINIMUM AV, THREE-IMPULSE MANEUVERS 


Run 

T r hr 

D, deg 

Plane changing 
associated with 
AVj, percent 

Energy changing 
associated with 
AVj, percent 

Plane changing 
associated with 
AVg, percent 

Energy changing 
associated with 
AV^, percent 

Plane changing 
associated with 
AV^, percent 

Energy changing 
associated with 
AV 3 , percent 

1 

8 

0 

-- 

58.4 

-- 

1.0 

— 

40.6 

2 

8 

2 

23.8 

67.3 

0.8 

. 1 

75.4 

32.6 

3 

8 

4 

69.2 

52.2 

4.4 

.6 

26.4 

47.2 

4 

8 

6 

71.4 

56.8 

14.5 

.9 

14. 1 

42.3 

5 

12 

0 

— 

63. 5 

— 

.8 

— 

35.7 

6 

12 

2 

57.4 

64.9 

14.0 

. 7 

28.6 

34.4 

7 

12 

4 

22.0 

65.3 

25.9 

.5 

52. 1 

34.2 

8 

12 

6 

43.4 

65.2 

15. 1 

15.6 

41. 5 

33.2 

9 

24 

0 


72.0 

-- 

. 5 

-- 

27.5 

10 

24 

2 

61.3 

72.2 

4.6 

.6 

34. 1 

27.2 

11 

24 

4 

49.4 

73.2 

7. 5 

.8 

43. 1 

26.0 

12 

24 

6 

64.4 

72.8 

6.8 

.7 

28.8 

26. 5 

13 

48 

0 

-- 

75. 3 

— 

.3 

— 

24.4 

14 

48 

2 

46.2 

76.9 

5.8 

.3 

48.0 

22.8 

15 

48 

4 

72. 1 

76.6 

6.6 

.4 

21.3 

23.0 

16 

48 

6 

64.2 

76.7 

12.0 

.4 

23.8 

22.9 

17 

72 

0 

— 

78.3 

-- 

.2 

-- 

21. 5 

18 

72 

2 

57. 1 

78.5 

7.8 

. 3 

35. 1 

21.2 

19 

72 

4 

67.3 

78.4 

10.3 

. 3 

22.4 

21. 3 

20 

| 72 

6 

57.4 

78.4 

17.4 

.4 

25.2 

21.2 


(Tj = 8 hours and D= 6°), AV does 32.1 percent of the plane changing and 
40. 8 percent of the energy changing; AVg does 56. 0 percent of the plane changing 
and 28.2 percent of the energy changing; and AVg does 11.9 percent of the plane 
changing and 31.0 percent of the energy changing. 

The AV components and generalized eccentric anomaly angles for the minimum 
AV, three-impulse maneuvers are given in table II in units of Earth radii and hours. 

It is observed that e ^ - 6 . — 180° and e g’ - e £ ^ 180° for all cases, with e ^ 

always small, positive, and acute. The flight time between first and third impulses 
for most maneuvers was within an hour of the stated period of the first transfer conic. 

As a further illustration, consider again computer run 4 (T^ = 8 hours and 

D = 6°) illustrated in figure 1. For this case, e ^ = 1.05°, = 199.78°, 

0 2 , = -160.50°, 0g' = -6. 17°, and eg" = -3.90°. The constraints for run 4 were well 

satisfied with converged values for the four components of the constraint vector and 

for the 11 components of the vector dF'/da of the order of 10 . Approximately this 

degree of numerical accuracy was attained in all cases examined. 
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TABLE II. - CONTROL-VECTOR AV COMPONENTS, GENERALIZED ECCENTRIC ANOMALY ANGLES, 
TOTAL AV, AND FLIGHT TIMES FOR THE MINIMUM AV, THREE-IMPULSE MANEUVERS 


Run 

Au 3 

Av 3 

Aw 3 

l h 

Au 2 

Av 2 

Aw 2 

0 X 

A“i 

4 V ! 

AWj 

F 

1 

0. 00000 

0.86246 

0. 00000 

5. 5850 

0.00000 

0.02083 

0. 00000 

5.6057 

0. 00000 

1. 2439 

0. 00000 

2. 1272 

2 

. 02204 

. 86894 

. 14545 

4. 7932 

-.00011 

.00467 

-.00478 

6.4110 

.02733 

1. 2401 

.20415 

2. 1451 

3 

.06396 

.84661 

. 18669 

4.3831 

.00250 

.02781 

-.05038 

6. 7632 

.03710 

1.2386 

. 24068 

2. 1892 

4 

. 11027 

.83462 

.21315 

4. 1818 

.01057 

.03707 

-.09788 

6.9075 

.04022 

1. 2378 

.25691 

2.2385 

5 

. 00000 

.73213 

. 00000 

6. 3988 

. 00000 

.01625 

. 00000 

6.4169 

. 00000 

1. 3759 

. 00000 

2. 1243 

6 

.01595 

. 72800 

.09653 

5.6231 

-.00047 

.01747 

-.01565 

7. 1844 

.01591 

1. 3736 

. 16148 

2. 1412 

7 

.04728 

. 72066 

. 13974 

5. 0302 

. 00264 

.02223 

-.04812 

7. 7397 

.02573 

1. 3720 

.20975 

2. 1769 

8 

.08001 

.71204 

. 164 70 

4.8726 

.00833 

.02892 

-.08644 

7. 8485 

.02713 

1. 3716 

.21894 

2.2160 

9 

. 00000 

. 58 103 

. 00000 

8.0704 

. 00000 

.01050 

. 00000 

8.0848 

. 00000 

1. 5294 

. 00000 

2. 1209 

10 

.01036 

. 57840 

.06403 

6.9593 

-.00011 

.01157 

-.01458 

9. 1860 

.00933 

1. 5279 

. 13060 

2. 1341 

11 

.02855 

. 57386 

.09425 

6.4962 

.00193 

.01440 

-.03881 

9.6176 

.01257 

1. 5273 

. 15331 

2. 1587 

12 

. 04716 

.56813 

. 11839 

6. 3827 

.00476 

.01833 

-.06451 

9.6955 

.01300 

1. 5272 

. 15708 

2. 1848 

13 

.00000 

. 488 10 

. 00000 

10. 175 

. 00000 

.00671 

. 00000 

10. 186 

. 00000 

1. 6240 

. 00000 

2. 1188 

14 

.00675 

.48631 

.04551 

8.8195 

. 00004 

.00745 

-.01189 

11.532 

. 00466 

1. 6232 

.09516 

2. 1285 

15 

.01745 

.48305 

.07113 

8. 4604 

.00109 

.00910 

-.02844 

11.868 

.00572 

1. 6230 

. 10598 

2. 1449 

16 

.02833 

.47859 

.09443 

8. 3720 

. 00240 

.01154 

-.04533 

11.931 

.00588 

1. 6230 

. 10769 

2. 1621 

17 

. 00000 

.45110 

. 00000 

11.650 

. 00000 

.00515 

. 00000 

11. 660 

. 00000 

1. 6617 

. 00000 

2. 1180 

18 

.00523 

.44962 

.03872 

10. 191 

. 00006 

.00571 

-.01008 

13. 110 

.00297 

1. 6612 

.07681 

2. 1259 

19 

.01314 

. 44677 

.06295 

9.8755 

.00073 

. 00694 

-.02307 

13. 407 

. 00352 

1. 6611 

.08402 

2. 1387 

20 

] 

.02118 

.44271 

. 08574 

9. 7969 

.00154 

.00883 

1 

-.03621 

13.465 

. 00362 

1. 6611 

.08514 

2. 1520 


CONCLUDING REMARKS 


It is believed that the numerical study makes several quantitative contributions. 

1. The study verifies numerically the existence of three-impulse maneuvers 
which afford measurable improvement in 18 of 20 cases studied, even for realistically 
small launch windows, over corresponding one-impulse maneuvers. For the cases 
examined, the improvement is from -0. 23 percent to +7. 40 percent lower velocity 
penalty, that is, total AV. 

2. For the cases examined, the study obtains the optimum relative distribution 
(among the impulses) of energy-changing and plane- changing duties. Previous studies 
have obtained minimum AV, three-impulse trajectories from sets of maneuvers in 
each of which all of the plane changing occurs at the high-altitude middle impulse. 
Imposition of this characteristic — which is valid for the minimum AV, infinite-time 
maneuver — would have been unrealistic for the (finite-time) cases examined, for 
which the low-altitude third impulse performed as much as 29. 3 percent of the plane 
changing. 
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3. The study determines that the velocity penalty for the small- launch- window 
three-impulse cases examined is only 0. 38 to 5.2 percent higher than for corresponding 
coplanar cases. 

4. The analysis discerns that the three- impulse maneuvers involve a 0. 13 to 

5. 8 percent greater velocity penalty than the corresponding absolute minimum AV, 
time-open three-impulse maneuvers. 

The program of the accelerated gradient method has a high degree of flexibility 
and can be adapted to a variety of parameter optimization problems. The current study 
was conducted on both IBM 7094 and Univac 1108 computers, and a typical case con- 
sumed 10 to 15 minutes of machine time. A large percentage of the core storage of the 
7094 was used. After the current problem was incorporated into the program, the only 
usual remaining difficulty was to select appropriate penalty constants for the penalty 
function phase of the convergence process and a first guess for an optimum control 
vector. A second computational phase inherent in the accelerated gradient program, 
using a Newton method for functional minimization subject to nonlinear constraints, 
encountered numerical convergence problems only infrequently. The accelerated 
gradient program is not usually strongly dependent on the first guess for the control 
vector in problems with a payoff function with reasonably uncomplicated contours. 
However, the level of complexity of the three-impulse study usually required correct 
guessing of at least the algebraic sign and order of magnitude of each component of the 
control vector before convergence would ensue. An occasional further difficulty was 
the obtaining of a degenerate solution, that is, one in which an impulse or coast dura- 
tion vanished, even though a nondegenerate solution existed. Degenerate solutions 
were rejected, and alteration of the first guess for the optimum control vector eventu- 
ally led to attainment of a true three- impulse maneuver. 

A further related study under consideration deals with minimum AV, two- 
impulse and three-impulse transfer onto a trans-Mars asymptotic velocity vector from 
a regressing oblate Earth-assembly parking ellipse, with fixed period of the first 
transfer conic. Another study underway seeks minimum AV, four-impulse transfer 
from an Earth parking ellipse onto a conjunction-class trans-Mars asymptotic velocity 
vector, with fixed time of flight between first and fourth impulses. 

Recommended for further study are (1) solution of specific fixed- time four- 
impulse problems, (2) construction of a good first-guess generator to use as an adjunct 
to the optimization method for fixed-time four-impulse problems, (3) determination of 
absolute minimum AV, three-impulse or four-impulse maneuvers for fixed-time prob- 
lems, and (4) modification of an optimization method to eliminate degeneracies of the 
type in which an impulse or a coast vanishes. 
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